Internal Josephson Oscillations for Distinct Momenta Bose-Einstein Condensates 
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The internal Josephson oscillations between an atomic Bose-Einstein condensate (BEC) and a 
molecular one are studied for atoms in a square optical lattice subjected to a staggered gauge field. 
The system is described by a Bose-Hubbard model with complex and anisotropic hopping parame- 
ters that are different for each species, i.e., atoms and molecules. When the flux per plaquette for 
each species is small, the system oscillates between two conventional zero-momentum condensates. 
However, there is a regime of parameters in which Josephson oscillations between a vortex-carrying 
atomic condensate (finite momentum BEC) and a conventional zero-momentum molecular conden- 
sate may be realized. The experimental observation of the oscillations between these qualitatively 
distinct BEC's is possible with state-of-the-art Ramsey interference techniques. 

PACS numbers: 67.85.Hj, 05.50.+q, 64.70.Tg, 67.60.Bc 
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I. INTRODUCTION 

Quantum coherence is a subject of fundamental im- 
portance and practical interest, especially concerning the 
construction of quantum logic devices. The Josephson 
effect, historically proposed to occur in superconduc- 
tors, has become an important tool for quantum coher- 
ence measurements. In cold atoms, a quantum super- 
position between two chemically different species (atoms 
and molecules), which is yet another ramification of the 
same effect, has been observed by means of Ramsey- like 
interference experiments for Bose-Einstein condensates 
(BECs) in a trap [H-Q. Later, experiments and theoret- 
ical studies included an optical lattice However, 
the studies were restricted to relatively simple quantum 
systems. 

The use of Berry's phase to realize artificial gauge fields 
for cold atoms has proven to be very fruitful to emulate 
more complex quantum systems 0- The possibility of 
creating vortex lattices with this technique @ provides 
an exciting prospect of reaching the fractional quantum 
Hall regime. Besides, the generation of a staggered mag- 
netic flux in a driven two-dimensional (2D) optical lattice 
holds promises of unprecedented simplicity [9j-fTo|. By 
loading a staggcred-flux optical lattice with cold bosonic 
atoms, distinct super/fluids can form, depending on the 
value of the flux per plaquette. For <f> < tt, the bosons 
condense at zero momentum, whereas for <j) > n a fi- 
nite momentum BEC is realized, which carries a vortex- 
antivortex lattice [ToL [Lsj . In addition, when manipulat- 
ing the interactions in the system by means of a Feshbach 
resonance, a bound state of two bosonic atoms (a Fesh- 
bach molecule) can occur, thus raising even further the 
parameter space for the realization of different BECs: in- 
deed, each type of particles, atoms (a = 1) and molecules 
(er = 2) can, in principle, condense either at zero or at 
finite momentum. 

Here, we first study the two-body problem of the 2D 
staggered- flux lattice for cold bosonic atoms. The stag- 



gered flux splits the lattice into A and B sublattices, thus 
introducing a pseudospin degree of freedom into the prob- 
lem. Due to the breaking of time-reversal symmetry and 
the presence of the pseudospin, bound states always ap- 
pear, irrespective of the repulsive interaction strength. 
This surprising finding provides a unique opportunity for 
realizing cold atoms experiments in which pseudospin de- 
grees of freedom play an important role in scattering pro- 
cesses. They may also shed light on the corresponding 
fermionic many-body problem in the context of high-T c 
cupratcs. Indeed, a staggered- 7r-flux phase was first pro- 
posed by Marston and Affleck [l7| to describe the pseu- 
dogap regime of high-T c cupratcs, and it has been ad- 
vocated by many to be the hidden-order behind high-T c 
superconductivity (l8| . 

Second, by taking into account a molecular formation, 
we use the Bogoliubov theory to study the collective be- 
havior of the generalized Bose-Hubbard model at zero 
temperature. The interplay between species and pseu- 
dospin degrees of freedom results in an effective four-band 
description that supports various out- of-p hase collective 
modes, also known as Leggett modes [19j . In particular, 
we find a regime of parameters in which coherent oscil- 
lations between two qualitatively distinct BECs can be 
realized: When the flux per plaquette for each species 
0cr=i.2 < I" and the hopping amplitudes J a =i,2 > 0, 
the system oscillates between two conventional zero- 
momentum condensates. However, for tt < 4> a —\ < 37r, 
3tt < (j) a =2 < 4-7T, and J a =2 < 0, a coherent oscillation 
between a conventional zero-momentum molecular BEC 
and a vortex-antivortex carrying atomic BEC can occur. 
The latter describes an internal Josephson oscillation be- 
tween two macroscopic groundstates carrying different 
quantum numbers [20( ■ 

The outline of this paper is the following: in Sec. II we 
introduce the model and calculate the two-atom scatter- 
ing, the many-body problem, and the collective modes 
in Subsecs. IIA, IIB, and IIC, respectively. In Sec. Ill 
we provide estimates for the realization of the collective 
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FIG. 1. (Color online) Optical square lattice in the presence 
of an artificial staggered flux, (a) Schematic representation of 
a staggered flux pattern with strength <j> in each plaquette. (b) 
Two interpenetrating square sublattices .4 and B connected 
by the vectors ei,l £ {1,2,3,4}. The unit cell is a square 
lattice (red circles) with a lattice constant a. 



modes. In Sec. IV we discuss the possibilities of ob- 
serving our results experimentally, and we present our 
conclusions. 



II. THE MODEL 

We consider a 2D Bose-Hubbard model for single- 
species cold bosonic atoms in the presence of a staggered 
flux, 

H = E (^- 1 ) , / 4 at6 r+ei +h.c.) 
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Here, J is the hopping amplitude between nearest- 
neighbor sites, (f> is the flux per plaquette, which al- 
ternates in sign between neighboring plaqucttcs, I € 
{1, 2, 3, 4}, and U is the on-site repulsive Hubbard inter- 
action. The operators aj(a r ) and bl(b T ) are the bosonic 
creation (annihilation) operators at site r in the square 
sublattices A and B, respectively, and n r is the num- 
ber operator. The sublattice constant is given by a = 
A/-\/2, where A is the laser wavelength and the connect- 
ing vectors e; are defined by ei = — e 3 = (A/2) e x and 
e2 = — e-4 = (A/2)e y ; see Fig. [TJ With the techniques 
described in Ref. a strongly driven time-dependent 
optical lattice can give rise to an effective hopping am- 
plitude J for the Hamiltonian ([1]), which takes positive 
as well as negative values [2l[. A physical estimation of 
the parameters based on Ref. [1] will be discussed in Sec. 
III. 



A. Two-atom scattering 



tical lattice. The two-body Schrodingcr equation related 
to the Hamiltonian (JTJ) is given by 



hx2+hx2^>H r2 + U6 rul2 



* = Et> (2) 



where 
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is the single-particle kinetic term for particle i = 1,2, 
with the discrete displacement operator A r in the pres- 
ence of staggered flux (/> defined as 



A r *(r) 



A^(r) = 



-j|e^/ 4 [*(r + ei ) + *(r + e 3 )] 
+e"^ /4 [*(r + e 2 ) + tt(r + e 4 )]|, 
- J je-^ /4 [*(r + ei ) + #(r + e 3 )] 
+e^ /4 [*(r + e 2 ) + *(r + e 4 )]l, (4) 



and I 2x2 is the 2x2 identity matrix. The two-particle 
wavefunction >F is constructed by taking the tensor prod- 
uct of two spinorial wavefunctions (>F_4(ri), "Fe( r i)) <8> 
A( r 2) , ^ s( r 2)) , or, in terms of its component, 

* T = [*AA(ri, r 2 ), *4e(ri, r 2 ), Vra(*i, r 2 ), * B s(ri, r 2 )] 
Similarly Eq. ([2]) can be written explicitly as 

/US rur 

V o 

By using a unitary matrix [22 
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Let us start by considering the scattering of two 
bosonic atoms at lattice sites ri and r 2 of the square op- 



Eq. ([5]) can be expressed in terms of the center-of-mass 
R = (ri + r 2 )/2 and relative r = iq — r 2 coordinates, 
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Using plane-wave states S^>(ri,r 2 ) = e' R e' r VK,k 
with center-of-mass and relative quasimomenta K and k, 
respectively, we rewrite the non-interacting Schrddingcr 
equation, which may then be diagonalizcd by using an- 
other unitary matrix S', to yield 



[E- A 



where A^- k is the i-th eigcncncrgy corresponding to the i- 
th component of the "pscudospin" eigenvector $o,K,k = 
S'^K,k- The subscript zero in the eigenvector ^o.K.k 
denotes the non-interacting limit of the problem. Since 
the interaction matrix (k|f/|k ) = Ul^i is diagonal and 
momentum-independent, the Lippman-Schwinger equa- 
tion for the scattering problem, 

|$K,k) = |$0,K,k) + G(£ K ,k)l>|$K,k) 

can be resummed to all orders in U in the new basis 
&o.K..k{ r ) for each pscudospin component to yield 
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*gk(r) 



= <K,k 
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UG^(E KM ,r) 
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Here, G^(E,r) is the Fourier transform of the non- 
interacting Green's function G«> (E, k') = l/(E-\\ Ck> + 
iO). In contrast to the usual atomic scattering problem, 
the pseudospin-carrying particles give rise to four inde- 
pendent scattering continua, sec Fig. [2J Each scattering 
continuum, in addition, may support a bound state ly- 
ing above it, which occurs when the scattering amplitude 
diverges, i.e., at the position of the pole with energy E 
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see solid curves in Fig. [2] 

We observe two striking features in the present model, 
which hold for all flux values except for the cases <f> = 
2im, n £ Z, when the two sublattices become degenerate. 
Firstly, the existence of at least one flat band A^ =0 k = 0, 
which is independent of the relative quasimomentum k 
at zero center-of-mass quasimomentum K = 0, guaran- 
tees the existence of a bound state solution with energy 
Eb = U, even for arbitrarily weak repulsive interaction 
strength. This is in stark contrast to the usual Hubbard 
model in a simple 3D cubic lattice, where a critical re- 
pulsive interaction strength is required for the existence 
of a bound state solution [23[ . Secondly, at certain ener- 
gies both scattering and on-resonancc processes can take 
place simultaneously in separate pscudospin channels. 
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FIG. 2. (Color online) Energy spectrum E (in units of 
the hopping parameter J) of the two-body Hamiltonian (TI]) 
as a function of the center-of-mass quasimomentum K for 
U/J — 1.818 and 4> = 11 ■ The four shaded bands I-IV depict 
the scattering continua for the four pseudospin components, 
and the solid curves are the energy dispersion of the repul- 
sively bound pairs that lie above their respective scattering 
continuum. 



Collisions between free atoms and repulsive bound pairs 
can thus lead to interesting dynamics since the pseu- 
dospin is generally not a good quantum number. More- 
over, due to the lack of dissipation in an optical lattice, 
the repulsive bound pairs are expected to be long-lived 
[23| . and it is conceivable that they can form a BEC. 



B. Many-body problem 

To study the dynamics when both atoms and compos- 
ite particles are present, we generalize the Hamiltonian 
((T|) to an effective two-species model: 

H = ~ E (^e i ^ ( - 1)!/4 4,a& r +e (l a+h.C.) 



E 

reA-:-B.a 



(e CT - /i CT )n r:Cr + ^TV^Tlr.cr - 1) 



^ a\ 2 ar,iarS + E K,2^,lb r ,l + h.C. 



,(H) 



with species index a = 1(2) denoting atoms (composite 
particles), and an inter-species conversion term charac- 
terized by the strength g. Here, e CT= i(2) and /v=i(2) are 
respectively the on-site energy and the chemical poten- 
tial for atoms (composite particles), and fi 2 = 2(J-i- The 



4 



hopping amplitude J ai the flux strength <p a , and the on- 
site interaction U a are species dependent because of the 
difference in masses m-i = 2m\. Multi-band effects are ne- 
glected because only small inter-atomic interactions are 
considered. In addition to repulsive bound pairs, compos- 
ite particles (which we henceforth call molecules) can also 
be formed by the use of Fcshbach resonances 24[ ■ In the 
absence of a staggered flux, this two-species model can 
in fact be derived from a microscopic Fcshbach resonance 
model in which the relative on-site energy e - =1 /e cr=2 an d 
the inter-species coupling strength g can be tuned via 
an external Feshbach magnetic field 0, [2fj|. The 
molecule-molecule interaction U2 and the atom-molecule 
interaction terms are neglected in the present work be- 
cause they are of higher order in the molecule operators, 
and we are considering the limit of a low molecule filling 
factor. Throughout our work, we consider ni <C 0.3 and 
these terms have a negligible contribution. 

The corresponding phase diagram exhibits a rich be- 
havior, which includes an array of superfluid-Mott- 
insulator transitions Q and a quantum Ising transition 
due to an enlarged symmetry group U{\) x Z 2 (27L l28l|. 
The Ising transition, which can also be viewed as a topo- 
logical confinemcnt-dcconfinement transition, separates 
a phase with a molecular condensate (MC) and nor- 
mal atoms with a residual Z2 symmetry, from a phase 
in which both molecules and atoms (MC+AC) are con- 
densed. The mean-field phase boundary can be deter- 
mined by generalizing the work of Ref. [5| and identify- 
ing the development of an instability in the atomic spec- 
trum in the phase with only a molecular condensate; see 
Fig. G3 We will henceforth measure the energy param- 
eters in units of ftio = 10kHz, which is the typical har- 
monic frequency of each well on a lattice site. At higher 
values of the interaction strength, the staggered flux is 
expected to simply rcnormalize the critical value (U/4J) C 
of the supcrfluid-Mott-insulator transition at a fixed den- 
sity, as it was found for single-species [l(J ■ 
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FIG. 3. (Color online) The phase boundary of the quantum 
Ising transition plotted as a function of the molecular con- 
densate fraction 712 and the coupling strength g/hu. Each 
lattice site represents a sufficiently deep harmonic trap with 
frequency kui = 10kHz, in terms of which all energy units are 
based. The figure is plotted for the following chosen param- 
eters: .h/hto = -0.11, J 2 /hLj = -0.08, Ui/huj = 0.2, ei/Hu = 
1.5, t2/hio = 2, 4>i = 4-7I-/3, and <f>2 = 167r/5. 




FIG. 4. (Color online) Collective modes of the atom-molecule 
system subjected to a staggered flux: (a) Scenario (I) with 
Ji/huo = 0.11, J2/huj = 0.08, and ^1,2 = 47r/5. (b) Scenario 
(II) with J-i/hw = -0.11,J 2 /ftu = -0.08, (j>i = 4n/3, and 
4>2 = I671-/5. The other parameters are chosen as follows: 



g/hu = 0.8, (7i /hw = 0.2, ei /hw 
and hu) = 10kHz. 



l.5,e 2 /Ml = 2, m = 0.5, 



Collective modes 



where 



Deep in the superfluid regime, however, BEC states 
with distinct momenta can concurrently form because 
the atomic and molecular fields are subjected to different 
flux values 4> a . In this subsection, we focus on the collec- 
tive modes in the superfluid regime where the U(l) x Z2 
symmetry is broken. By performing a canonical trans- 
formation to the band representation, 



k, cr 



£k.o 



1 

y/2 V l^k.rr 

1 / £k, c 
V^V |ek,c 



(12) 



fflk.cr + &k,cr I 



£k,, 



4J 



4 J V 2 



kx't 



a \ . / k v ia 
sin I — — sin 1 



(13) 



is the kinetic form factor in the momentum space (the 
prime denotes sublattice coordinates), we find that the 
non-interacting part of the Hamiltonian (|11[) is described 
by four distinct bands with energy dispersions 



±2jJ cos 2 [fc+a] + cos 2 [AT 



-2 cos[(/) (J /2] cos[fc + a] cos[fc a] 



1/2 



(14) 
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where fc* = (k x / ± k y ')/2 (see Fig. 2 for the definition 
of the x' — y' coordinate system). The on-site interac- 
tion and atom-molecule coupling terms take the general 
inter-/intra-band coupling form. We then consider the 
formation of BECs, where the lowest energy operators for 
the atoms (a = 1) and for the molecules (a = 2), both 
acquire a non-zero expectation value at the respective 
condensation points ko(a = 1,2), with condensate num- 
ber iVo,<T. The Bogoliubov approximation amounts to 
substituting a~ =kaia) a ^ ^/N^ + a k=ko{a)a and keep- 
ing the fluctuations up to the quadratic order. The 8x8 



grand-canonical mean-field Hamiltonian then reads 



iuj n ,k 

where the matrices 



<5k = g 
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M fci2 = 
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if/ini(u_ fc )* + g^/n^v^ k 

\Uinxiul)* +gy/n^v k 




\U 1 n 1 u k + g^/n^{4Y 


±U 1 n 1 u k h +g v / n^{v k )* 
( £ k,iY, 



-fc,2 



'fc,2^ ' 



with the Nambu spinor defined as 



*fc,i ' 



k,l' a fc+fco(o-=l),2' a ^fc+fc ((T=l),2' a fc+fco(cr=l),2' 



fc+fc ((T = l),2 



i^a; n + + e a - \i a + 2U a n a , where u„ is 

4 = 1± '7fe,l } 7-fc,l»?fco(<r=l),l J 7fco(<T=l),l. w fc 



Here, ni(n 2 ) denotes the atomic (molecular) condensate fraction, 
the Matsubara frequency. The transformation coefficients are given as: 

1±? 7o,2 ? 7fc.i V-k.i, s k = l±V*k +ko (a=i)a 7 lko(cr=i),iVk,i, where r}k, a = Arg{e^" /4 cos[(A; x < + k y >)a/2] + e ~^" /4 cos[{k x > 
fcj,/)a/2]} and ko(a) denotes the quasimomentum at which species a condenses. 

I 



We require that the terms that are linear in the fluc- 
tuations vanish, yielding the relations, 

sign(j 2 ) exp(?4 ((7=2))2 + 2?7 Mct=1):1 ) = 1, 
sign(j 2 ) exp(77* o((T=2)i2 + ri ko{a = 1)A 

+ %o(o-=2)-fe (o-=l),l) = 1 ' ( 16 ) 

as well as 

1 + sign( J 2 ) exp(r?fc o(o , =2)i2 + 2t? m<7=1)>1 ) 

XgfM V / ^2 4o(cr=2),2fco(a-=l) 

+n CT (l + <5 CT , 2 ) (^ CT - /v + e* + U a n a ) = 0, (17) 

where 5 a .b is a Kronecker delta. We find that the two- 
species system is only allowed to have two condensa- 
tion scenarios: (I) for J 2 > 0, ^ 0i >2 < tt and (II) for 
J 2 < 0, 7r < 0i < 37r, 37T < 2 < 47r. In scenario (I), both 
atoms and molecules condense at the same momentum 
ko = 0, whereas in scenario (II) the molecules condense 
at k = 0, while the atoms condense at a different mo- 
mentum k = (ir/a,Tr/a). Outside these regimes, the 
system does not admit a self-consistent stable mean-field 



condensate state. Since /x 2 = 2/xi, we can equate the two 
expressions resulting from solving Eq. (|17[) for fi a and 
obtain a relation between n± and n 2 , demonstrating that 
the two condensate fractions are not independent. 

The excitation spectrum for both scenarios, obtained 
by solving the mean- field Hamiltonian (|T5j) . consists of 
four collective modes because of the A— B sublatticcs 
and the atom-molecule species degrees of freedom; see 
Fig. (U]). The lowest branch is the Goldstone mode asso- 
ciated with both, A— B sublattices and atom- molecule in- 
phase density fluctuations. The next lowest branch cor- 
responds to the A—B out-of-phase but atom- molecule in- 
phase oscillation mode. The third and fourth branches, 
separated by a gap from the two lowest branches, de- 
scribe collective modes where atoms and molecules oscil- 
late out-of-phase. In fact, the long-wavelength physics of 
the third branch is equivalent to a coherent oscillation of 
density number of the two species in real time. It thus 
describes an internal Josephson effect of an atomic and 
a molecular condensates that can have distinct quantum 
numbers, with a frequency that is given by the energy 
gap. The dependence of the gap on the coupling strength 
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FIG. 5. (color online) Josephson frequency of two BECs with 
distinct momenta as a function of coupling strength for dif- 
ferent on-site interaction strength, with hw — 10kHz. The 
system exhibits dynamical instabilities for g < 0.7ftw. 



g for scenario (II) is shown in Fig. (|5|) for weak and in- 
termediate on-site interactions. This result provides a 
promising experimental opportunity to study the coher- 
ent dynamics of two distinct macroscopic ground states 
in a cold atomic system. 



III. EXPERIMENTAL OBSERVATION 

The many-body phenomena discussed in Sec. II. C may 
be observed experimentally by using the set-up proposed 
in Refs. 0, [21(, where ultracold bosons arc trapped in 
a time-dependent square optical lattice with staggered 
currents. Here, we provide an estimate of the parameters 
specific for this configuration. Using the Floquet theory, 
one can show that the time-dependent problem can be 



described by an effective Hamiltonian of the form [21 



■I'irh-'Mlrl 



-JT+-i-K[M,N]+H int 



(18) 



where f + = YlreA,l{ a * br + e i + h - c -} and 
[M,N] = 2J2 re4l {(-l) l alb r+ei - h.c.} are the 
single-particle hopping operators, with I = 1,2,3, 
and 4. Jo(z) and J\(z) are Bcsscl functions of the 
first kind and H m t is the onsitc interaction term. 
The parameters J = J dxdy w(x,y)w(x + 1/2, y), 
£at = 2kVq J dxdy \w(x, y)\ 2 cos(kx) cos(fcy), and 
£m = kVq J dxdyw*(x + A/4, y)[sin 2 (fca;) — 
cos 2 (ky)]w(x — X/4,y) arc given in terms of the 
Wannicr function w(x,y) of the optical lattice with 
amplitude k and lattice depth Vo. The time-periodic 
driving has a frequency f2, and the strongly driven 
regime is fulfilled by the conditions J,£,m,U <C Ml. 
For simplicity, we may take the Wannier function to be 
the ground state wavefunction of a harmonic oscillator 



w(x,y) = (l/l H oV^)cxp[-(x l/^HO\_ 
site with harmonic length Iho = (£'r/Vb) 1 ^ 4 (A/27r), 
where E r = h 2 k 2 /2m is the recoil energy, m is the 
mass of the atom, and k = 2tt/X. For an optical lattice 
potential with amplitude fixed at k = 1 and Vo/E r ps 2, 
we get ^jy — 1-404 Vq. For the other parameters, we 
use the result from a full band calculation made in Ref. 
[2l| that gives J ~ 2.5 ^m- Because the Bessel functions 
range from negative to positive values, we see that the 
tunneling amplitudes J and K can attain negative values 
as a function of the driving frequency, see Fig. 6(a). 

In the superfluid regime of the Hamiltonian (fT8|) , there 
exist four distinct order parameters in four ranges of driv- 
ing frequency, which we denote here as superfluid A, B, 
C and D; see Fig. 6(b): superfluids A and C are the zero- 
momentum condensate while superfluids B and D are the 
finite- momentum condensates. As explained in Ref. [2l| . 
even though two of the four superfluids carry the same 
momentum, they have different symmetries associated 
with them. On the other hand, there is a square region 
enclosing the origin of the phase diagram that character- 
izes a Mott insulator phase, not shown in Fig. 6(b). For 
the lattice amplitude we consider here, the width of the 
Mott region at unit filling is ~ 0.04J7, thus not sizable, 
and the system does not cross into the Mott insulator 
state as the driving frequency is varied. 

To relate to the Hamiltonian studied in Eq. (1), we 
may rewrite the kinetic term of Eq. (|18p in the polar 
form 



■E 

reA,l 



(Jeffe^-l)'/"^ 



h.c. 



(19) 



with the effective parameters given as 



</eff 



Wn 2 



g6v 
Ml 



26v 
Ml 



off = 4 tan" 



J Jo ' 



(20) 



In this representation, the negativity of the tunneling 
amplitude arises from the phase value e m ' . 

An additional ingredient of the model we study in this 
work is that it consists of two species, namely atoms 
and molecules, trapped in the same time-dependent op- 
tical lattice. We take the simplest consideration where 
their masses are related by m2 = 2m\. Accordingly, the 
parameters of the model, through their dependence on 
the Wannier functions, scale as Iho,<j=i = 1-189 Iho,<j=2, 
J a =i = 3.838 J a =2, €n,o-=i = 0.901 ^N,a=2, and J a=2 ~ 
2.5£m.o- = 2- The additional subscript cr, as previously 
used, denotes the species index. Thus, for the molecules 
subjected to the same time-dependent optical lattice, the 
positions of the zeros of the tunneling amplitudes J a =2 
and K a= 2 (the real and imaginary components, respec- 
tively) are both scaled by a factor 0.901 as compared to 
the atomic case in Fig. 6(a). Consequently, the driv- 
ing frequency range for the various molecular superfluid 
phases is shifted accordingly. 
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FIG. 6. (color online) (a) The value of J/ J (the real com- 
ponent) and K/J (the imaginary component) plotted as a 
function of 2^n/M1, for amplitude k = 1 and Vo/Er = 2. 
The four zones A, B, C and D characterize four superflu- 
ids with different symmetries in the weak interaction regime, 
(b) Phase diagram of the Hamiltonian ()18p . The solid curve 
shows the superfluid phases that the system realizes as the 
parameter 2£jv/ftfi is tuned from to 8.5. We take U/J ~ 1 
for our case. 



In Sec. II. C we propose two possible scenarios to study 
coherent oscillations of atomic and molecular conden- 
sates. In the phase diagram of Fig. 6(b), the two scenar- 
ios are: (I) both atoms and molecules condense at zero 
momentum and form the superfluid A; (II) the atoms 
condense at finite momentum (superfluid B), whereas 
the molecules condense at zero momentum (superfluid 
C). Therefore, for the setup under consideration in this 
section, we estimate that for a driving frequency range 
< 2£jv.<r=i/ft^ < 1, scenario (I) can be realized, 
whereas for the driving frequency region 2£jv CT =i/^ ~ 3, 
scenario (II) can be realized. 



IV. DISCUSSIONS AND CONCLUSIONS 

Bound states with pseudospin components can be ob- 
served experimentally, since the binding energy can 
be inferred by spectroscopic measurements [IE He EmI ■ 
The lifetime of the bound state can also be probed 
through the techniques used in Ref. [HI • More impor- 
tantly, a coherent oscillation between distinct macro- 
scopic ground states can be experimentally observed by 
performing double-pulse Ramsey experiments [l|, H| for 
bosonic atoms in the optical lattice setup described in 
Ref. The regime of negative hopping parameters 

J < for the staggered flux optical lattice may be real- 
ized experimentally, as shown in Ref. |3lj . Recent ex- 
periments in optical lattices without flux have probed 
the negative hopping regime by shaking the lattice with 
a periodic perturbation |32j , similar to the one occurring 
in Ref. @. 

In conclusion, we investigated the dynamics of a gen- 
eralized two-species Hubbard model subjected to a stag- 
gered flux. At the two-atom level, we studied the scatter- 
ing problem with an onsite interaction. The local nature 
of the interaction allowed us to solve the two-body prob- 
lem exactly, and we obtained the scattering spectra, as 
well as bound states induced by the repulsive interaction. 
The fact that the atoms carry a pseudospin gives rise to 
new scattering properties that are not present in conven- 
tional cold atomic systems. In particular, we find that 
a repulsive bound state always exists at zero center-of- 
mass quasimomentum due to the existence of a flat band. 
We then use the existence of bound states in the two- 
body problem as a motivation to study the many-body 
problem in a generalized Hubbard model subjected to a 
staggered flux. The resulting system can execute various 
collective density oscillations, especially an out-of-phase 
mode, which is also known as Josephson oscillation. Fur- 
thermore, we propose a way to realize a quantum super- 
position of different species with different momenta, zero 
for the molecular condensate and finite for the atomic 
one. The results can be probed by combining state-of- 
the-art techniques. This work opens up new perspectives 
in the realization of more complex BEC's and in the pre- 
cise control of their dynamics. 



ACKNOWLEDGMENTS 

We acknowledge financial support from the Nether- 
lands Organization for Scientific Research (NWO). We 
are grateful to A. Hemmerich and I. Spielman for a crit- 
ical reading of this manuscript. 



[1] E. A. Donley, N. R. Claussen, S. T. Thompson, and C. [3] R. A. Duine and H. T. C. Stoof, Phys. Rep. 396, 115 

E. Wieman, Nature (London) 417, 529 (2002). (2004). 

[2] S. J. J. M. F. Kokkelmans and M. J. Holland, Phys. Rev. [4] N. Syassen et al, Phys. Rev. Lett. 99, 033201 (2007). 

Lett. 89, 180401 (2002). [5] D. B. M. Dickerscheid et al, Phys. Rev. A 71, 043604 



(2005). [19 
[6] M. J. Bhaseen, A. O. Silver, M. Hohenadler, and B. D. [20 

Simons, Phys. Rev. Lett. 103, 265302 (2009). [21 
[7] J. Dalibard, F. Gerbier, G. Juzeliunas, and P. Ohberg, [22 

arXiv:1008.5378 (to be appeared in Rev. Mod. Phys. Col- 

loq.). [23 
[8] Y.-J. Lin, R. L. Compton, K. J. Garcia, J. V. Porto, and [24 

I. B. Spielman, Nature 462, 628 (2009). 
[9] A. Hemmerich and C. Morais Smith, Phys. Rev. Lett. [25 

99, 113002 (2007). 
[10] Lih-King Lim, C. Morais Smith, and A. Hemmerich, [26 

Phys. Rev. Lett. 100, 130402 (2008). [27 
[11] Lih-King Lim, A. Hemmerich, and C. Morais Smith, 

Phys. Rev. A, 81, 023404 (2010). 
[12] N. R. Cooper, Phys. Rev. Lett. 106, 175301 (2011). [28 
[13] N. R. Cooper and Jean Dalibard, EPL 95, 66004 (2011). 
[14] E. Kapit and E. Mueller, Phys. Rev. A 83, 033625 (2011). [29 
[15] M. Aidelsburger, M. Atala, S. Nascimbene, S. Trotzky, 

Y.-A. Chen, and I. Bloch, eprint arXiv:1110.5314. [30 
[16] G. Moller and N. R. Cooper, Phys. Rev. A 82, 063625 

(2010). [31 
[17] I. Affleck and J. B. Marston, Phys. Rev. B 37, 3774 

(1988). [32 
[18] S. Chakravarty, R. B. Laughlin, D. K. Morr, and C. 

Nayak, Phys. Rev. B 63, 094503 (2001). 



A. J. Leggett, Prog. Theo. Phys., 36, 901 (1966). 

A. J. Leggett, Rev. Mod. Phys., 73, 307 (2001). 

A. Hemmerich, Phys. Rev. A, 81, 063626 (2010). 

J. Sabio, F. Sols, and F. Guinea, Phys. Rev. B 81, 045428 

(2010). 

K. Winkler et al., Nature 441, 853 (2006). 

C. Chin, R. Grimm, P. Julienne and E. Tiesinga, Rev. 

Mod. Phys. 82, 1225 (2010). 

R. B. Diener and T.-L. Ho, Phys. Rev. Lett. 96, 010402 
(2006). 

H. P. Buchler, Phys. Rev. Lett. 104, 090402 (2010). 
L. Radzihovsky, J. Park, and P. B. Weichman, Phys. Rev. 
Lett. 92, 160402 (2004); Leo Radzihovsky, Peter B. We- 
ichman and Jae I. Park, Ann. Phys. 323, 2376 (2008). 
M. W. J. Romans, R. A. Duine, S. Sachdev, and H. T. 
C. Stoof, Phys. Rev. Lett. 93, 020405 (2004). 
C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Na- 
ture (London) 424, 47 (2003). 

H. Moritz, T. Stoferle, K. Giinter, M. Kohl, and T. 
Esslinger, Phys. Rev. Lett. 94, 210401 (2009). 
A. Eckardt, C. Weiss, and M. Holthaus, Phys. Rev. Lett. 
95, 260404 (2005). 

A. Zenesini, H. Lignier, D. Ciampini, O. Morsch, and E. 
Arimondo, Phys. Rev. Lett. 102, 100403 (2009). 



